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We discuss the alpha-muon sticking coefficient in the muon-catalysed "in flight" d-t fusion in 
the framework of the Constrained Molecular Dynamics model. Especially the influence of muonic 
chaotic dynamics on the sticking coefficient is brought into focus. The chaotic motion of the muon 
affects not only the fusion cross section but also the fi — a sticking coefficient. Chaotic systems lead 
to larger enhancements with respect to regular systems because of the reduction of the tunneling 
region. Moreover they give smaller sticking probabilities than those of regular events. By utilizing a 
characteristic of the chaotic dynamics one can avoid losing the muon in the /iCF cycle. We propose 
the application of the so-called "microwave ionization of a Rydberg atom" to the present case which 
could lead to the enhancement of the reactivation process by using X-rays. 
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I. INTRODUCTION 

The muon catalyzed fusion (/iCF) of hydrogen iso- 
topes, especially d-t fusion, has been studied as a realiz- 
able candidate of an energy source at thermal energies. 
In the liquid D2 and T2 mixture, the muon assists the 
fusion through the formation of a muonic molecule, since 
the size of the muonic molecule is much smaller than that 
of the ordinary molecules and the fusing nuclei tend to 
stay closer. After the fusion process the muon is released 
normally and again it is utilized for another fusion. The 
efficiency of the /zCF is governed by the muon-sticking 
on the a particle which is produced in the fusion • 
The muon is lost from the /iCF cycle due to sticking (wo), 
unless it is not released successively through the inter- 
action with the medium. The rate of the stripping of 
the stuck muon from the a particle is known as the re- 
activation coefficient R and thus the effective sticking 
probability^;^) is determined by 
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wo(l - R). 



(1) 



The determination of the value of R is discussed in the 
reference 0. In this paper we do not take into account 
the medium effects which are supposed to be important 
to determine the precise value of R in the actual exper- 
imental setup. We rather aim to propose a method in 
order to enhance the reactivation process, by making use 
of the stochastic instability of the stuck muon in an os- 
cillating field. For this purpose we are mainly interested 
in investigating the impact of the regular and chaotic 
dynamics [4| . The experimental value of the initial stick- 
ing loq which is determined assuming the theoretical R 
in is tabulated in the reference 0. In the table the 
values of ujq from 7 separate measurements are smaller 
than theoretical estimate for the most part. The direct 
measurement of ujq has been conducted as well and gave 



the w = 0.69 ± 0.40 ± 0.14 % 0. The temperature de- 
pendent phenomena in the muon cycling rate and in the 
muon loss probability, which is a function of ufj^ , have 
been reported on the ddfj, at the temperature from 85K 
to 790K and on the dtfi from 5K to 16K 0,0 by mea- 
suring the fusion neutron yield and the K a X-ray yield 
lately. They have observed that the muon loss probability 
increases and the muon cycling rate decreases as the tem- 
perature decreases. In the latter case u e Jf varies from 
0.64% to 0.86±0.01% as the temperature varies from 16K 
to 5K. The temperature dependence in the muon loss 
probability is thought to be caused by a change of the 
reactivation coefficient R in Eq. Q 7, 8]. 

At thermal energies, where the /zCF takes place, fluc- 
tuations might play a role. We investigate the influence 
of the fluctuations by using a semi-classical method, the 
constrained molecular dynamics (CoMD) approach 0. 
As it is well known the molecular dynamics contains all 
possible correlations and fluctuations due to the initial 
conditions(events). In the CoMD, the constraints restrict 
the phase space configuration of the muon to fulfill the 
Heisenberg uncertainty principle. The results are given 
as an average and a variance over ensembles of the quan- 
tity of interest which is determined in the simulation. 
Especially we determine the enhancement factor of the 
reaction cross section by the muon as a function of the 
incident energy in the "in flight" fusion [To| : 
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where Q =17.59 MeV is the decay Q-value of this re- 
action. The enhancement factor of each event indicates 
the regularity of the system. The reaction is known 
to give a smaller reaction rate with respect to the fusion 
through the formation of the molecular complex: 

tn+D 2 -> [(tdfi)dee\* -> 4 He ++ +/i+n+d+2e+Q. (3) 

(for review dlH). 

Subsequently we determine the initial muon sticking 
probability, using the phase space distribution of the 
muon at the internal classical turning point. The stick- 
ing probability is evaluated regardless of the difference 
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of the two reactions eq. J5J) and 10 due to the fact that 
the Q-value is much larger than the binding energies of 
muonic atoms or molecules. We suppose that loq is insen- 
sitive to the formation process of the muonic molecules. 
A distinctive feature of our study is that we do not as- 
sume the ground state of the muonic molecule as the 
initial state of the muonic molecules, instead we use the 
initial state configuration by simulating the fusion pro- 
cess employing the imaginary time method 0, Q, 0] . 
As a consequence, in fact, the muon does not stick nec- 
essarily to the ground state of the alpha particle and 
this fact plays an important role when we proceed to 
the stripping of the bound muon in the oscillating field. 
The chaotic dynamics could prevent the muon from being 
lost in the /iCF cycle due to the sticking. It is achieved 
by utilizing the characteristic as a nonlinear oscillator of 
the trapped muon on the alpha particle. We draw an 
analogy between the muonic He ion in the present case 
and microwave-ionization of Rydberg atoms Il7l Il8| , 
where the driven electron in the highly excited hydro- 
gen atom in a strong microwave electric field exhibits the 
chaotic dynamics and is ionized. Since highly excited 
states in the atom, with high quantum principle number 
n, are in the quasi-classical regime, its stability can be 
explained in classical mechanics in terms of resonances. 
We carry out a numerical simulation by enforcing an os- 
cillating field (linearly polarized, oscillatory electric field) 
on the system. This can be, likely, achieved by radia- 
tion of a coherent Synchrotron Orbital Radiation(SOR) 
X-ray experimentally. The oscillating force causes the 
resonance between the force itself and the oscillating mo- 
tion of the muon around the alpha, especially when the 
driving frequency coincides with integer multiples of the 
eigen frequency of the muonic helium. In other words 
the muon can be stripped by controlling the chaos of the 
system. In passing we mention that an attempt to apply 
magnetic fields for the purpose of stripping from excited 
states of muonic helium has been proposed in |l9j |. 

This paper is organized as follows. In Sec. [H] we de- 
scribe the theoretical framework of the CoMD for muonic 
molecule formation and following fusion process briefly. 
The relation between the enhancement factor and the 
chaotic motion of the muon is discussed in Sec. IIIII We 
develop in Sec. IIVI a formula to estimate the initial a-/i 
sticking probability (u>o) and determine loq. Sec. 0is de- 
voted to the discussion of a possibility of muon release. In 
Sec. IVII we summarize the paper and mention the future 
perspectives of this study. 



II. FRAMEWORK 

The details of the framework of the CoMD is discussed 
in the references 0, |2(], ■ In the following we sketch 
the framework briefly by applying it to the case of the 
reaction J5J. We assume the ground-state t[i as targets at 
the beginning of the collision. The ground-state muonic 
tritium configuration in the phase space is obtained us- 



ing the CoMD approach |21| . Denoting the position of 
the particles (i = [i,t) in the phase space by (ri,pi) and 
the relative distance and momentum between fj, and t 
by r^t = |r M — r t | and p^t = |p M - Pt| respectively, the 
modified Hamilton equations for the muonic tritium with 
constraints are 
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(4) 
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where we use relativistic kinematics; £i = y^pfc 2 + mfc 4 
and U(r^) — —e 2 /r fJit is the potential of the muon. The 
second term in the eqs. (@J and © represent the con- 
straints: the phase space density of two distinguishable 
particles should be always equal to 1 or less |23, i.e., 
the terms prevent the muonic tritium from collapsing. 
In terms of the phase space distance, this can be writ- 
ten as r^ t p^ t > where £h = 1 and the equal sign 
should be satisfied in the case of the ground state muonic 
atom. For this purpose A^, the Lagrange multiplier for 
the Heisenberg principle, is determined depending on the 
distance of the muon from the triton in the phase space 
r^tP^t- If r^tP^t is (smaller)larger than £ H h, Af is a pos- 
itive (negative) small finite number. The approach gives 
the average binding energy of the ground state muonic 
tritium atom BE tf2 = —2.73 keV. The value is in agree- 
ment with BE H (M ll /M e ) = -2.71 keV, where BE H is 
the binding energy of the hydrogen atom, and M e and 
Af M are the reduced masses of the the electronic and 
muonic atoms. 

Using the obtained ground-state configuration as an 
initial state, we perform the numerical simulation of the 
fusion process (0 using 



Pi = -V r Z7(rj) 



(6) 



for all the particles(z = t,d,fi). As the interaction 
we consider a modified Coulomb interaction U(r) — 
Sj(^) **/ r y x (! ^ e ~ bTti )i with qt and q 3 being the 
charges of the particles and b — 9500A _1 . 

In order to treat the tunneling process in the frame- 
work of the molecular dynamics, we define the collective 
coordinates R co " and the collective momentum p coU as 
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oil — 
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Pd - Pt, 



(7) 



with T t ,Td (pt,Pd) being the coordinates (momenta) of 
the triton and the deuteron, respectively. To obtain the 
classical turning points, we first simulate the elastic col- 
lision. Subsequently we repeat the simulation in the tun- 
neling region by switching on the collective force, which 



is determined by F 



coll 



ycoll 



and 



-pcoll 



ycoll 



to 



enter into imaginary time [Til Il5l| . We follow the time 
evolution in the tunneling region using the equations, 
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dr 
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dr 



= -V r U(r? id) )-2FZ\, (8) 
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where r is used for imaginary time to be distinguished 
from real time, and p^ rfl are position and momen- 
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turn of the triton (the deuteron) during the tunneling 
process respectively. 

We assume that the fusion process occurs at small im- 
pact parameters compared with the radius of the muonic 
atom, i.e., we carry out the simulation of the head on 
collisions. Under this assumption, the penetrability of 

the barrier is given by miia 



11(E) = {l + exp(2A(E)/h))' 
denoting the action integral A(E) as 



A(E) 



(9) 



(10) 



with r a and r\> being the classical turning points. The in- 
ternal classical turning point is determined using the 
sum of the radii of the target and projectile nuclei. Sim- 
ilarly from the simulation without muon, we obtain the 
penetrability of the bare Coulomb barrier H (E). We 
choose the initial inter-nuclear separation 3 A. This is 
much larger than the scale of the muonic tritium radius, 
which is of the order of 1.3 mA. 
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FIG. 1: Enhancement factor by the bound muon (top panel) 
and A/^// M (bottom panel) as functions of the incident 
center-of-mass energy. The arrows in the figure indicate the 
point where total energy is zero. 



III. ENHANCEMENT OF THE CROSS 
SECTION BY THE MUONIC SCREENING 
EFFECT 

We introduce the enhancement factor of the cross sec- 
tion by the bound muon / M 

/ju = <T(E)/a (E), (11) 

where a(E) and <tq(E) are the screened cross section 
and the bare cross section, respectively. We approximate 
Eq. (jllfl by taking the ratio of the penetrabilities in the 
presence and in the absence of the muon: 

f fl =U(E)/U (E). (12) 

In the following discussion, the enhancement factor is 
referred as an indicator of the regularity of the muonic 
motion HH3. It plays a role of a sort of order param- 
eter and it is determined through the values obtained in 
the numerical simulation. In the top panel of Fig. ^ 
we plot / M as a function of the incident center-of-mass 
energy between the triton and the deuteron. From our 
simulation of the collisions using an ensemble of events, 
we determine the average enhancement factor and its 

variance: A/ M = \f* — fp 2 ] 1 ^ 2 . These are shown by 
squares and error-bars, respectively. Both the average 
fp and its variance increase exponentially as the incident 
energy decreases. The dashed and dotted lines in the 
figure correspond to the enhancement factor / M in the 
conventional and exact adiabatic limit respectively (See 
Appendix) . The average is in good agreement with the 



exact adiabatic limit with a screening potential = 8.3 
keV. This fact suggests that the enhancement mainly de- 
pends on the energies involved, i.e., the energy difference 
of the initial and the final atoms. Thus the problem 
of fusion is somewhat independent on the detail of dy- 
namics. Even if we had prepared the entrance channel 
through the formation of the molecular complex eq. J2J , 
the actual fusion cross section would not differ from the 
one calculated here. In particular in the limit E <C 
the fusion cross section can be written from eq. 111(1 and 
ea. (|A.2|l as 

Hid ME) = Hm ^{E + UJ - ^M e ^ u ^ , (13) 

-C/ — Hj hi — >U Up 

where r](Up) is the Sommerfeld parameter. 

In the bottom panel the ratio Af 2 /fp versus incident 
energy is plotted. In the high energy limit the ratio 
approaches zero, i.e., the /„ distribution becomes a 5- 
function (A/ M = 0) and the average / M approaches 1: 
there is no effective enhancement. In the low energy limit 
Af 2 /fp 3> 1, which implies that the system exhibits a 
sensitive dependence of the dynamics on initial condi- 
tions, i.e., occurrence of chaos. It is noteworthy that the 
slope of the ratio A/^// M changes at the ionization en- 
ergy of the muonic tritium, which we indicated by the 
arrows in the figure. At this incident energy the total 
energy of the system is zero. The total system might be 
unbound at the incident energies higher than this point, 
while the 3-body system is bound at lower energies. We 
indeed verify the manifestation of chaos by plotting the 
Poincare surface of section with respect to the enhance- 
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ment factor for two events in Fig. [2] In the figure we 
show the surface of section for two selected events at the 
incident energy 0.18 keV on the x-p x plane (FIG. left 
panels) and on the z-p z plane (FIG. |3 right panels), re- 
spectively. We choose the beam axis to coincide with 
the z-axis. At the incident energy 0.18 keV the average 
enhancement factor, / M = 1.1 xlO 27 as one can see in 
Fig.[T] In the top panels, with / p = 4.1xl0 19 (< f^) and 
the ratio of the external classical turning point in the 
presence of the muon(rj) to the one in the absence of 
the muon(r°): r£/r a = 0.15, the points show a map of a 
typical regular event. By contrast in the bottom panels, 
with / M = 2.7xl0 31 (> ffj,) and the ratio of the classical 
turning points r^/r° a = 0.06, the points show the map of 
an irregular event: the points cover a large section of the 
map. The irregular muonic motion leads to smaller ex- 
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FIG. 2: Surface of section for 2 events, one has small 
(top panels) and the other has large / M (bottom panels), on 
the x-px(leR panels) and the z-p z (right panels) planes at the 
incident c.o.m energy 0.18keV, in the atomic unit 

ternal classical turning point. As a consequence it gives 
larger enhancements factors opposed to theprevious re- 
sults in the case of the electron screening 0, |2(J , where 
the irregular(chaotic) events give smaller enhancement 
factors. This contradiction is accounted for the fact that 
the system is bound in the present case at low incident 
energies, while in the previous case even the lowest inci- 
dent energy which was investigated is much higher than 
the binding energy of the electrons. Therefore the chaotic 
dynamics of the electrons causes to dissipate the kinetic 
energy between the target and the projectile and lowers 
the probability of fusion. 

We can easily deduce the reaction rate(A dt ) at the liq- 
uid hydrogen density(pLH = 4.25 x 10 22 cto~ 3 ) and at 
very low thermal energies fcgT <C U^, 



A 



dt 



= lim a{E)\i t - r d \p LH (14) 
= lim a (E)f ll \i t - r d \p LH . (15) 



Using the fact that the average enhancement by the muon 
is written as = ao(E + U^/a^E) in terms of the 
screening potential U, L and substituting the data in the 
NACRE compilation 23] for the bare reaction cross sec- 



tion CT (t/ M 



eq.JT3J), 
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<//, 



2U U 



M, 



dt 



PLH ~2.8x lOV 1 ], 



(16) 



where we denote the reduced mass of d and t as Mdt and 
the charge symmetry in the t+d system has been taken 
into account as it is explained in the appendix. This value 
is in agreement with the result obtained with an indepen- 
dent approach [l(J], X dt — 0.5 x lQ 5 [s -1 ]. We, neverthe- 
less, stress that this value is obtained using the average 
enhancement factor. As we have seen in Fig. El fluctua- 
tions might be extremely large and there are events which 
have the / M more than 10 4 times larger than the average 
f^. Therefore the "in flight" d—t fusion rate could fluctu- 
ate as well and become comparable to the inverse lifetime 
of the muon. We will discuss the role of fluctuations more 
in detail in a future work. 



IV. MUON STICKING PROBABILITY 

We estimate the sticking probability of muons on the 
alpha particle in the exit channel: 



(td/j,y 



'He + M ' 

/4, 



( 4 He/ir 



(17) 
(18) 



Recall that the Q-value is 17.59 MeV, while the ener- 
gies involved both in the "in-flight" reaction JSJl or with 
molecular complex formation eq. @ is of the order of sev- 
eral eV, thus we expect that the actual entrance chan- 
nel for fusion is irrelevant to the following dynamics of 
tunneling, fusion, decay and sticking of the muon to the 
alpha particle. 

The muon remains bound (eq. H18fl h if the binding 
energy of the muon on an alpha particle: 



BE a u = 



M a , 



2e 2 



(19) 



is negative, in the center-of-mass system of the muon and 
the alpha particle. We denote the reduced mass of p. and 
a as M afl . The effect of the finite nuclear mass must be 
taken into account, because a muon is 206.8 times heavier 
than an electron. From this condition, BE af ^ < 0, we 
deduce the following equation for the angle 9 between 
and r n . 



M, 



2c- 



COS0 > 



^(|r,| 2 + |r a | 2 ) 
M QM |r M ||r Q | 



(20) 



The condition Eq. (|20|l is fulfilled when g, the r.h.s of 

the equation, is equal to 1 or less and for the solid angle 
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n = 2tt(1 —g) [steradian] in the 3-dimensional space. We 
can therefore estimate the sticking probability by H/4tt, 
if g < 1. We point out that |ro,| in the equation is written 
MS i\ function of the decay Q-value: 



r„ = 



■y/2Q/M a 



(21) 



where r an = r n — r a is the relative velocity between 
the a and the neutron and M an is their reduced mass. 
In particular, the sticking probability can be estimated 
easily in some limiting cases: 
CASE 1. if IrJ is about |rj, 



9 



1 - 



2e 2 



|i>|/M ajU < 1 



(22) 



Furthermore, in addition, we assume that the muon is 
bound in the ground state of the 5 He at the moment of 
the fusion(adiabatic limit), i.e., — , 2 f r | = -10.942x2.0 

keV, g is estimated to be 0.89 and thus the sticking prob- 
ability is 5.6 %. In passing we mention that g is 0.352 
for the reaction d + d ^ 3 He +n + 3.268 MeV under the 
same assumptions with CASE 1. We deduce 32.4 % of 
the muon sticking probability in this case. In the case of 
the reaction t + t — > a + n + n + 11.33 MeV 24], if we 
assume that two neutrons bring away the maximum en- 
ergy 9.44 MeV, we can estimate g — 0.84 and the sticking 
probability= 7.84 % for this reaction. 
CASE 2. if r,, I is much smaller than |r„|, 



iJ*2l>l 
2 IrJ " ' 



(23) 



where we assume 



0.0. This means there is 



no sticking probability in this case. 

We remind that muons can have higher velocity com- 
ponents in the quantum mechanical system. Turning to 
the case of our simulations, the above velocity for the 
bound muon is obtained as the average velocity over the 
ensemble of events. By inspecting each event, one could 
find some events which satisfy the condition g < 1. In- 
deed one of the two events shown in the top panels in 
Fig. |3 which is regular, has g = 0.93(< 1.0) there- 
fore the sticking probability of this event itself is not 
zero (3.8 %). While the other in the bottom panels has 
g = 1.9(> 1.0), the sticking probability is zero. In the 
same way we calculated g for all the events which are 
created in our simulation. The resulting sticking proba- 
bility of the muon on the a particle is shown with filled 
circles in the top panel of Fig. |3| as a function of the 
incident energy of the collision. At the same time, we 
carry out the simulation of the exit channel by creat- 
ing 20000 events (randomly chosen directions of outgo- 
ing particles). The time integration of the equation of 
motion, Eq. © for all the particles(i = a,n,/j,) is per- 
formed using Predictor-corrector integration scheme. We 
distinguish the muon sticking event (|18fl from the release 
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FIG. 3: Incident energy dependence of the sticking probability 
of the muon on the a particle. The statistical error is shown 
by error bars, otherwise it is within the size of the points in 
the figure, (top panel) Distance between the muon and the 
alpha particle as a function of the inter-nuclear separation 
(bottom panel) 



event l|17[) by monitoring the binding energy of the muon 
on alpha particles Eq. I|19l) and the radius of the muonic 
ion, TZ at _, = |r M — r Q |. We count the events where BE a ^ 
and JZ afi maintain to be negative and small respectively, 
up to the point where the alpha is distant enough from 
the neutron. 

In the bottom panel in Fig. [3] typical trajectories of 
the distance between the muon and the alpha particle 
are shown as a function of the inter-nuclear separation. 
Among the 6 curves shown in the figure, 3 curves, which 
show oscillational behaviors, correspond to the sticking 
events. While 3 other curves increase monotonically af- 
ter R na exceeds 0.01 A. The horizontal straight line in 
the figure indicates the size of the radius of the ground 
state muonic He atom. The obtained sticking probabil- 
ity is shown with open circles with error-bars in the top 
panel in Fig. [3] First, as one can see clearly, the result of 
numerical simulation agrees with the sticking probability 
which is calculated considering the solid angle il. In the 
figure we plot the stickingcoefficient obtained from the 
direct measurement of luq |2| with a solid line with an er- 
ror bar. The resulting sticking probability range nearby 
the experimental value as a function of the incident en- 
ergy except for several points, which have zero and rel- 
atively large sticking probability. This fact supports our 
assumption that coq is insensitive to the formation pro- 
cess of the muonic molecules. We mention especially that 
the incident energies at which the sticking probability be- 
comes zero is slightly below than the ionization energy of 
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the muonic tritium. 



V. MUON STRIPPING 

Before concluding the paper, we would like to suggest 
that the stuck muon is possibly stripped from the alpha 
particle, by enforcing a linearly polarized oscillatory elec- 
tric wave on the system. The periodic motion of the 
stuck muon can be expressed in terms of nonlinear oscil- 
lations. For a nonlinear oscillator the oscillating driving 
force, i.e., linearly polarized, in the direction z, oscillat- 
ing field of frequency 7 and a peak amplitude F of the 
field, 



corresponds to 0.1 lnm and 0.44nm, respectively, in terms 
of the wave length. Such a X-ray is available from Syn- 
chrotron Orbital Radiation(SOR) facilities. In our nu- 
merical simulation, instead of the discontinuous frequen- 
cies, we get the proper frequency of the stuck muon for 
each event. The frequency is obtained by Fourier trans- 
form of the oscillation and with thus obtained frequency 
7^ we perform the simulation of stripping with above 
external force. In Fig. 0] we show the time T map of 
the oscillational motion of the muon on the r — p plane, 
where r is relative distance between the muon and the al- 
pha and p is its conjugate momentum. The time T map 
is obtained by sampling the points in the phase space at 
discrete times Hal 



zF sin(7t), 



(24) 



t = nT (T = 2tt/7, n = 1, 2, ■ • ■)• 



(25) 



causes the resonance between the force itself and the os- 
cillating motion of the muon at driving frequencies which 
are integer multiples of the fundamental frequencies of 
the muon. The same concept has been applied in the ion- 
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FIG. 4: Time T map of a stuck muon on the r-p plane(left 
panels) . One with the external oscillating force (right panels) 
with frequency 7, both in the atomic unit 

ization of the Rydberg atoms in a microwave field. There 
the highly excited atom, Rydberg atom, is prepared by 
laser excitation. In the present case of muons, the muonic 
He is not stuck necessarily in its ground state. So that 
the muon can be ionized directly from one of such an 
excited state by the radiation of an electrostatic wave, 
otherwise the muon is, at first, prompted to an excited 
state and then ionized. In either cases it will be achieved 
using X-rays, since the fundamental frequencies for the 
ground and the first excited state of the muonic He ion 



We follow 1000 cycles of the driven oscillation in our sim- 
ulation. We choose two sticking event which are shown in 
the bottom panel in the figure^ one has a smaller ampli- 
tude, a tightly bound state(top panels), and the other has 
a larger amplitude, a loosely bound state(bottom panels). 
The left panels show the map of the stuck muon without 
external force. The map remains in a limited manifold 
around r — 0. The right panels show the case with an 
external force with driven frequency 7 = 2 x 7^, where 
7 M = BE afl /h is the angular frequency of the muon in the 
muonic helium ion, i.e., we are investigating the 2:1 reso- 
nance. In the case of the tightly bound muon, the muonic 
atom is excited in a loosely bound state and then ionized. 
While the tightly bound muon is ionized in one step. The 
muonic atom in the external oscillating field is captured 
into the resonance and ionized due to its stochastic insta- 
bility [2(| . One can see clearly that the muon is expelled 
from the Helium with the external oscillational force with 
the corresponding frequency of the unperturbed system. 
We point out that the /mM molecule is not destroyed by 
the external force with same frequency which we used in 
the above discussion. 



VI. CONCLUSIONS AND FUTURE 
PERSPECTIVES 

In this paper we discussed the alpha-muon sticking co- 
efficient in muon-catalysed "in flight" d-t fusion. We per- 
formed numerical simulation by the Constrained Molecu- 
lar Dynamics model. Especially the influence of muonic 
chaotic dynamics on the sticking coefficient is brought 
into focus. The chaotic motion of the muon affects not 
only the fusion cross section but also the /1 — a sticking 
coefficient. The irregular (chaotic) dynamics of the bound 
muon leads to larger enhancements with respect to reg- 
ular systems because of the reduction of the tunneling 
region. Moreover they give smaller sticking probabilities 
than those of regular events. 
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We proposed a method to strip the stuck muon from 
the alpha particle by exposing the system in the X-ray 
radiation field. Its numerical experiments have been per- 
formed under an oscillating external force with the driv- 
ing frequency twice as high as the angular frequency of 
the stuck muon and the muon has been released success- 
fully with the selected frequency. By utilizing the chaotic 
dynamics one can prevent the muon from being lost in 
the /LtCF cycle due to sticking. 

Based on these results, in our future study, we will de- 
velop a theory to investigate the temperature dependent 
phenomena, including the sticking probability and the 
muon cycling rate, which are reported by experimental- 
ists 0, 0, 0] • Further quantitative analysis of the muon 
stripping with the oscillating force should be undertaken. 
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APPENDIX: ENHANCEMENT FACTOR IN THE 
ADIABATIC LIMIT 



In the section lllll we introduced the enhancement fac- 
tor of the cross section by bound muon in terms of the 
barrier penetrability. In the case of the electron screen- 
ing, one often assumes that the effect of the screening can 



be represented by a constant shift, 
barrier and replaces ea. fTTl) by 



Uu, of the potential 



lit 



ME) ' 



(A.l) 



V ' y, is called screening potential. (Jq(E) can be rewritten 
in terms of the S-factor S(E) and Sommerfeld parameter 
7](E), by writing down the incident energy dependence 



of the barrier crossing rates explicitly [27J , 

hi 

In the limit of <C E, is approximated by 
fu = ex P 



n V (E)^ 



(A.2) 



(A.3) 



The merit of this conventional formula is that one can 
easily estimate the upper limit of the enhancement by 
using the adiabatic approximation in the framework of 
the Born-Oppenheimer approximation. In the present 
case of the muonic tritium target, reflecting the charge 
symmetry in the t+d, is obtained by assuming equally 
weighted linear combination of the lowest-energy i "ger- 
ade" and "ungerade" wave function for the muon. 



fu 



2 [ 



a Trr](E) - 



■)• 



(A.4) 



it is given by substituting U u = BE t — BEs^ ~ 8.3 
keV with BE t and BE^ He being the binding energy of 
the ground state muonic tritium and the ground state 
muonic 5 He respectively and = 39.6 eV. The point 
is, however, that the procedure is justified only in the 
limit where U u is much smaller than E. is much 
larger than the low incident energies of our interest. The 
dotted curve in Fig. ^ which corresponds to Eq. (|A.4|I , 
indeed overestimates the enhancement with respect to 
the exact formula: 



fu 



(A.5) 



where we approximated nothing but S(E + U e ) ~ S(E). 
To apply this formula in the present case again one needs 
to take a linear combination of the muonic molecular 
states. This is taken into account and gives the dashed 
curve in Fig. ^ The curve is in accord with the average 
enhancement factor obtained from our simulation. 
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